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Abstract 

We investigated the formation and evolution of CDM subhalos in galaxy-sized and cluster-sized halos 
by means of ./V-body simulations. Our aim is to make clear what the "dwarf galaxy problem" is. It has 
been argued that the number of subhalos in simulated galaxy-sized halos is too large compared to the 
observed number of dwarfs in the local group, while that in cluster-sized halos is consistent with observed 
number of galaxies in clusters such as the Virgo cluster. We simulated nine halos with several different 
mass resolutions and physical scales. We found that the dependence of the cumulative number of subhalos 
iV c on their maximum circular velocity V c is given by N c oc V~ , down to the reliability limit, independent 
of the mass of the main halo. This implies that simulations for cluster-sized halos give too many halos with 
V c ~ 140km/s or less. Previous comparisons of cluster-sized halos gave much smaller number of subhalos in 
this regime simply because of their limited resolution. Our result implies that any theory which attempts 
to resolve the missing dwarf problem should also explain the discrepancy of the simulation and observation 
in cluster-sized halos. 

Key words: cosmology: dark matter — cosmology: theory — galaxies:clusters:general — galaxies: dwarf 
— mcthods:n-body simulations 



1. Introduction 

The Cold Dark Matter(CDM) scenario (White & 
Rees 1978) has been the standard theory of the formation 
and evolution of the structures in the Universe. In this sce- 
nario, galaxies and clusters of galaxies are formed bottom- 
up. It has been remarkably successful in explaining the 
large scale structures (Davis et al. 1985) and numerous 
observational results (e.g. Springel et al. 2005). Until re- 
cently, however, the small-scale structures of CDM, like 
subhalos in a galaxy-sized halo, could not be studied by 
numerical simulation because of the lack of the computa- 
tional power. Recent improvement of computational pow- 
ers made it possible to study not only dark matter halos, 
but subhalos in a parent halo. Usually these subhalos are 
interpreted as corresponding to galaxies in cluster-sized 
halos and satellite dwarf galaxies in galaxy-sized halos. 

Klypin et al. (1999b) and Moore et al. (1999, here- 
after M99) reported that about 1000 subhalos formed in a 
simulated galaxy-sized halo. The number distribution of 
subhalos as the function of the circular velocity normal- 
ized by that of the parent halo turned out to be remark- 
ably similar for galaxy-sized and cluster-sized halos. This 
result, of course, is the direct outcome of the scale-free 
nature of the gravity and almost power-law amplitude of 
the density fluctuation. 



This similarity, which looks quite natural from the the- 
oretical side, showed serious discrepancy with observa- 
tion. The number of subhalos in the simulated CDM halo 
agreed well with the observation of the cluster of galaxies, 
while the simulation of a galaxy-sized halo predicted too 
many dwarf-sized halos. This discrepancy is now known 
as the "dwarf galaxy problem" , or "missing-dwarf prob- 
lem" , and has been one of the main topics in the study of 
galaxy formation and structure formation in the Universe. 

Solutions so far proposed for this problem can be clas- 
sified into the following two categories: 

• Dark matter is different from the CDM model in 
small scales. 

• Only a small fraction of subhalos is observed as 
dwarf galaxies. 

Many modified dark matter models had been proposed. 
Self Interacting Dark Matter(Spergel & Steinhardt 2000) 
is a CDM with small self-interaction cross section, 
which dumps small-scale fluctuation. Warm Dark 
Matter (Bardeen et al. 1986) has a power spectrum which 
has cut off at small scale. Some simulations with 
these dark matter models have been done. (e.g. Dave 
et al. 2001, Colin et al. 2000) 

If we assume that the standard CDM is correct and 
consider the possibility that not all subhalos are observed 
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as dwarf galaxies, the question is which subhalos are ob- 
served as dwarf galaxies. Stoehr et al. (2002) argued that 
observed dwarf galaxies are presently most massive sub- 
halos. by comparing numerical result with observation of 
velocity dispersion profile of dwarf galaxies. On the other 
hand, Kravtsov et al. (2004, hereafter K04) tracked each 
subhalo's evolution and showed many low-mass subhalos 
were massive in the past. Additionally, they introduced 
the model of star formation condition, and claimed that 
the number and distribution of satellite dSph galaxies was 
in agreement with that of such previously massive subha- 
los. 

Before jumping to such a solution, however, we should 
understand what is really the problem. In previous stud- 
ies, in order to determine the number of subhalos with 
circular velocity ~ lOkm/s, the detection limit was set to 
a few tens of particles. It is not clear whether or not we 
can rely on the result obtained with such a small num- 
ber of particles. Both softening parameters and the two- 
body relaxation might significantly affect mass and den- 
sity profiles of subhalos (e.g. Moore et al. 1996, Diemand 
et al. 2004b). Ghigna et al. (2000) showed that the low 
mass end of the circular velocity distribution of subhalos is 
affected by numerical resolution by comparing the results 
of two cluster-scale simulations with different resolutions. 
Same tendency can be seen in the more recent studies 
(e.g. Diemand et al. 2004a, Gao et al. 2004 and Reed 
et al. 2005). This means that the agreement in a cluster 
scale described in M99 might be a numerical artifact. 

In the present paper, as a first step to understand the 
formation and evolution of subhalos, we investigated the 
reliability of size distribution of subhalos in simulation 
results. We performed large-scale A-body simulations of 
formation of CDM halos. We simulated both galaxy-sized 
and cluster-sized halos. 

We carried out A-body simulations of nine halos of dif- 
ferent masses and resolutions and analyzed the distribu- 
tion of subhalos. From the results of these simulations, 
we estimated the effect of the mass resolution on the dis- 
tribution of subhalos. When the number of particles in 
one subhalo is small, both the mass and the circular ve- 
locity of that halo are affected, resulting in the deviation 
from simple power-law dependence. We formulated the 
reliability criterion for mass and circular velocity. 

We then compared the velocity distribution function ob- 
tained by A-body simulations with that of observed galax- 
ies in the Virgo Cluster. We found that the observational 
result is in good agreement with our medium-resolution 
result, which is affected by small- A effect around V c ~ 
lOOkm/s. In other words, our high-resolution simulation 
predicts the overabundance of moderate-sized galaxies. 

Previously the mass function of galaxies in a cluster 
has been believed to be well reproduced by CDM simula- 
tions (e.g. M99). Our result shows that the agreement, 
in the range of V c = 140km/s, is due to the lack of mass 
resolution in previous simulations. The discrepancy be- 
tween the simulated halos and observed galaxies extends 
to higher mass than in the galaxy scale. This is simply 
because there was no satellite of that mass in galaxy-scale 



calculations. Our result indicates that an analogue to the 
"missing dwarf problem" exists also in clusters of galaxies. 
Therefore, the "solutions" for the missing dwarf problem 
should give a distribution consistent with observations in 
both the galaxy scale and the cluster scale. 

The contents of this paper is as follows. In section 2, 
we describe the models and method used in our simula- 
tions. We give detailed explanation of the subhalo detec- 
tion method we developed. In section 3, we discuss the 
reliability limit of the number count of subhalos. Finally 
we summarize our result and discuss possible solutions to 
the missing dwarf problem in section 4. 

2. Model and Method 

2.1. Cosmological A -Body simulation 

2.1.1. Initial condition 

We adopted the Standard CDM model(SCDM), so 
that we can directly compare our result with that of 
M99 who adopted SCDM. Cosmological parameters are 
O m = 1-0, = 0.0, H a = 50km/s, ct 8 = 0.7. To generate 
the initial condition, we used the standard re-simulation 
method (Navarro et al. 1996). First we simulated 15 co- 
moving Mpc radius using 1.1 x 10 6 particlcs(m = 8.9 x 
10 8 M Q ) and 150 comoving Mpc radius using 1.1 x 10 6 par- 
ticles (m = 8.9 x 10 11 M Q ) to select candidates for galac- 
tic and cluster halos from the distribution of particles at 
z = 0. To select halo candidates, we used the standard 
Friends of Friends method (FoF, Davis et al. 1985) with 
the linking length which is 0.2 times the mean interparti- 
cle distance. Then we selected particles inside the radius 
four times larger than radius enclosing particles detected 
by FoF as the region of high resolution calculation, and 
traced these particles back to the initial condition. We 
replaced these particles with high resolution particles and 
re-ran the simulation. In this way, the external tidal force 
from outside the high resolution region was correctly taken 
into account. 

To generate initial conditions, we used the grafic2 pack- 
age (Bertschinger 2001). The calculation was done on 
an IBM pSerics690 of the Data Reservoir Project 1 of the 
University of Tokyo, and on a workstation with AMD 
Opteron 242 processors and 16GB memory. 

We prepared initial conditions of one region with three 
different resolutions for both galaxy-scale and cluster-scale 
runs. These runs are named Cl-H, Cl-M, and Cl-L and 
Gl-H, Gl-M, Gl-L, for cluster scales and galaxy scales. 
Here, postfixes H, L and M denotes high, medium and low- 
resolution runs, respectively. We also selected two regions 
for additional cluster-scale runs with medium resolution. 
These runs are named C2-M and C3-M. The initial red- 
shift z st art and particle mass m are listed in table 1. 

2.1.2. Time integration method 

The time integration was done in physical time and 
physical coordinates. In high-z regime, we varied time 
step At and the softening length e to reduce the time 

http:/ /data-reservoir. adm.s. u-tokyo.ac.jp/ 
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Table 1. Parameters for calculations. 



run 


m(M n ) 


At /?; „(10 b yr) e/i„(kpc) z start 




Gl-L 


4.1 x 10 b 


4.3 




0.6 63.1 


11.8 


Gl-M 


1.2 x 10 b 


2.9 




0.4 69.4 


13.1 


Gl-H 


5.1 x 10^ 


2.2 




0.3 74.0 


14.0 


C1-L,C2-L 


4.1 x 10 y 


6.5 




6 31.8 


5.56 


C1-M,C2-M 


1.2 x 10 y 


4.3 




4 36.7 


6.54 


C1-H,C2-H 


5.1 x 10 s 


3.2 




3 40.3 


7.27 


Table 2. Properties of halos. 


run 


M 


200 (Af ) 


A 20 o 


i? 2 oo(kpc) 


V c (km/s) 


N su b 


Gl-L 


2.3 x lO 1 ^ 


570490 


344 


236 


220 


Gl-M 


2.2 x lO 1 ^ 


1768684 


334 


235 


601 


Gl-H 


2.1 x lO 1 ^ 


4179435 


334 


233 


1346 


Cl-L 


1.7 x 10 lb 


405729 


3067 


1703 


350 


Cl-M 


1.4 x 10 lb 


1198496 


2934 


1495 


977 


Cl-H 


1.6 x 10 ib 


3170198 


3044 


1669 


2019 


C2-L 


8.2 x 10 i4 


199993 


2423 


1397 


129 


C2-M 


8. 


1 x 10 i4 


673493 


2416 


1373 


428 


C2-H 


8.2 x 10 i4 


1600368 


2424 


1203 


1010 



integration error. We adopted the procedure similar to 
that used in Kawai et al. (2004). We used the following 
formulation 



At(z) 
e(z) 



%crit H~ 1 

z+l 

%crit ~\~ 1 

z + l 



At 



fini 



-e fi 



(1) 

(2) 



where z cr n is the redshift at which we switch from z- 
dcpendcnt timcstep and softening to constant timcstcp 
and softening. The constant values are Atfi n and e/i n . 
These equations imply At oc \J\jp and e oc a for z > z CT n , 
where p is the mean density, a is the scale factor. The ini- 
tial redshift (z star t) and z cr u correspond to the standard 
deviation of density fluctuation (p — p)/p is 0.2 and 1.0, 
respectively. We chose Atfi n so that the energy conser- 
vation is better than 3% when calculations are done only 
for the high resolution particles. We chose e/i n so that it 
is sufficiently small compared to the size of subhalos. For 
runs with different N, we varied e/i„ in proportion to JVs, 
These values are summarized in table 1. 

For actual time integration, we used a parallel im- 
plementation of Barnes-Hut tree algorithm (Barnes & 
Hut 1986; Makino 2004) on GRAPE-6 (Makino et al. 2003) 
and GRAPE-6A (Fukushige et al. 2005). The opening an- 
gle was 0.3 for all run. Time integration was done using 
the standard leapfrog. 

Table 2 lists virial radii (-R200), masses (Af2oo) and num- 
bers of the particles (A 2 oo) inside i?200i and circular 
velocities (V c ) of the main halo. Here, -R200 is the radius in 
which the average density is 200 times the critical density. 
Figure 1 shows snapshots of a galaxy and a cluster halos 
from high resolution runs. 



2.2. Method for subhalo detection 

To detect subhalos in simulated galaxy halos, we 
developed a subhalo detection program based on the 
Hierarchical Friends of Friends (HFoF) method (Klypin 
et al. 1999a). The HFoF method is an extension of 
the standard FoF. The FoF algorithm with the linking 
length h detects the regions with the density higher than 
Ph = 2m/(^nh 3 ). In the HFoF method, FoF is repeat- 
edly applied to the regions found by FoF, with decreasing 
linking length, so that we can identify all density peaks. 
Initially, at level 1, we use the linking length 



hi 



1 2m 



2p c 



(3) 



Here pcrit is the critical density of the Universe and de- 
fined as pern = 3Ho/8irG. 

For halos detected at level i — 1, we apply the FoF 
method at level i with linking length 

h . =2 -(i-W hl} ( 4 ) 

i.e., we shrink the linking length by 2~ 1 / 3 at each level. 
This factor must be small enough so that the result does 
not depend on its value. When FoF of level i is applied 
to a halo found in level i — 1, there are three possibili- 
ties. The first one is that a single halo is again found, 
but with a smaller number of particles. Second is that 
nothing is found, or there are too few particles in the de- 
tected halo. The third possibility is that multiple halos 
arc found. When one or more halos arc found, we con- 
tinue to level i + When the halo vanished at level i, we 
regard that a single halo is detected at level i — The 
lower limit for the number of particles in a subhalo is 10. 
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Any method for subhalo detection based on the iso- 
density contour has a problem that many groups found 
are not gravitationally bound. To select gravitation- 
ally bound groups, we employed the evaporative method 
(Pfitzer et al. 1997). Using this method, we can discard 
unbound groups. In addition, we can remove unbound 
particles from bound groups. Since the iterative calcula- 
tion of gravitational potential necessary for this method 
is computationally expensive, we used GRAPE-6 to cal- 
culate the potential energy in order to accelerate calcula- 
tion. We applied this method to all subhalo candidates, 
even when the number of particles is very large. Thus, for 
all halos, we can accurately determine the bound mass. 



To summarize, our halo finding algorithm is follow- 
ing. First, the HFoF method is applied to all particles to 
construct density hierarchy. Then we apply evaporative 
method to each of detected density peaks (at the highest 
level for which that halo is found). If the peak is rejected 
by the evaporative method, that peak is regarded as not 
detected and we apply the same procedure to its parent 
halo (halo in one level lower). If the parent halo of one 
halo contain no other halo, we also regard that it is de- 
tected only at level i — 1. We apply this procedure to all 
hierarchies of subhalos, and regard all remained halos as 
detected. 

Total numbers of detected subhalos N su b are listed in 
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tabic 2. Right panels of figure 1 show the detected sub- 
halos. 

3. Results 

3.1. Cumulative distribution of subhalos 

Top two panels of figure 2 show the cumulative number 
N C (V C ) of subhalos as the function of the circular velocity 
V c for Gl and CI halos. In both cases, N c is proportional 
to V~ 3 for large subhalos (higher V c ), but for small sub- 
halos the number of subhalos drops off the power-law line. 
If wc investigate the curves carefully, we can see that N c 
drops off from the power law in two stages. For example, 
the curve for Gl-H starts to deviate from the power-law 
line at V c ~ 20km/s, and then show the second leveling off 
at V c ~ 6km/s. All other curves show similar two-stage 
behavior. 

To see this tendency more clearly, in bottom panels of 
figure 2 we plot the cumulative number normalized by 
V~ 3 . In Gl runs, runs H and M show narrow flat region, 
where the relation N c oc V~ 3 holds. The L (low-resolution) 
run shows hardly any flat region. Similar tendency is vis- 
ible for CI runs, though somewhat less clear. 

In these panels, it is clear not only the second leveling- 
offs but also the first deviations are dependent on the res- 
olution. In other words, first deviations are numerical 
artifact, and simulation results are reliable only for the 
subhalos larger than these first bendings. 

This result indicates that the region where the numer- 
ical result is reliable is quite narrow. For Gl-H run, the 
numerical result is reliable only for V c > 15km/s, for which 
only ~ 100 halos exist. 

Figure 3 shows the same cumulative number of subha- 
los, but now as the function of bound mass of halos M. 
We can see that the result is quite similar to that of figure 
2, except that the second leveling-off is not visible. 

The reason why the leveling-off is visible for velocity and 
not for mass is simple. Since we have posed minimum 
number of particles for subhalos (N > 10), there is no 
subhalo with N < 9, which is the reason why the curve 
terminate suddenly for N C (M) plot. In the case of N C (V C ), 
halos with TV > 10 but with small V c were counted, but 
halos with N < 10 are not counted no matter how high its 
V c is. Thus, N c levels off gradually. 

The bottom two panels of figure 3 show the cumulative 
number normalized by M _1 . Again, the bending is clearly 
visible. The "flat" region is not very clear, perhaps simply 
due to small number statistics. For example, if we inter- 
pret the result of run Cl-M naively, we might conclude 
that normalized plot in the mass range 10 11 to 10 12 is not 
flat. However, since there are only 20 subhalos with mass 
larger than 10 12 , this result is not statistically significant. 

For the small-mass region where runs with different 
mass resolutions show different results, the statistical 
noise is small because the number of halos is large. At 
the first sight the result in this region might look reliable, 
simply because of this low noise. However, clear difference 
between runs with different resolutions indicate that re- 
sult in this region is numerical artifact and is not reliable 



at all. 

From figure 3 we can estimate the minimum number of 
particles in the subhalos above which the number count is 
reliable. For run Cl-H, the mass at which the normalized 
cumulative count bends off from the flat line is between 
5 x 10 10 Mq and 10 u Mq, which correspond to the number 
of particles in the subhalos between 100 and 200. For 
other five runs, bending points are all similar, between 100 
and 200. So we can conclude that the cumulative count is 
only reliable for subhalos with number of particles larger 
than 200. 

In order to see this tendency more clearly, in figure 4 
we show the cumulative number of subhalos plotted as 
functions of the number of particles in subhalos. Thin 
solid curves are fitting function of the form 



N c (n p ) = 




where C and no are fitting parameters. For all curves, we 
use no = 220. The other parameter C was chosen by eye. 
We can see that for all runs formula (5) gives excellent fit. 

Reed et al. (2005) performed simulations similar to 
what we performed, and measured the cumulative num- 
ber of subhalos as functions of circular velocity and mass. 
They have not reported this bending at n p — 200. For 
velocity distribution, they did not see this bending sim- 
ply because their cutoff number of particles for subhalos is 
much bigger than what we used. As the result, the second 
leveling off occurs almost simultaneously as the first bend- 
ing. For mass distribution function, their figure 11 shows 
the same bending behavior as we see in figure 4, though 
the tendency is less clear simply because of the way the 
data are plotted, and they did not discuss the shape of 
the mass distribution function. Qualitatively same ten- 
dencies were reported by Gao et al. (2004) in galaxy scale 
and Diemand et al. (2004a) in cluster scale. 

3.2. True number of subhalos and comparison with the 
observation 

For our largest calculation with 4M particles in the 
virial radius of the parent halo, we found that the dis- 
tribution of subhalos is reliable only up to 100-200 most 
massive subhalos. For smaller subhalos, the dependence 
of the cumulative number on both the circular velocity 
and mass becomes significantly shallower than the sim- 
ple power-law. Thus, we should assume that the number 
of subhalos formed in real CDM cosmology would show 
this simple power-law behavior down to mass or circular 
velocity much lower than the limit for which the current 
simulations can give reliable result. 

The classic work by M99 used around 1M particles for 
the parent halo, and compared the distribution of around 
1,000 subhalos. On the other hand, we can conclude that 
their result is actually reliable only for around 20 most 
massive halos, and the number of smaller subhalos was 
significantly underestimated. 

Figure 5 shows the cumulative number of subhalos for 
runs C2-H, M and L, as well as that estimated for the 



G 



KASE AND MAKINO 





Fig. 2. Top: Cumulative number of subhalos as a function of V c . Bottom: Cumulative number of V c normalized by V c ~ 3 . Left and 
right are for Gl halos and CI halos. Solid, dashed, and dotted curves indicate the results of high, medium, and low resolution runs, 
respectively. 



Virgo cluster in the way same as in M99. We can see that 
run C2-M, which used the number of particles similar to 
that in M99, shows the best agreement with the obser- 
vation. The result of run C2-H shows the slope steeper 
than the observation starting at V c = 200km/s. The ac- 
tual number count for run C2-H becomes higher than ob- 
servation only at V c = 140km/s or around, because the 
number of massive halos in run C2-H is slightly smaller 
than the observational value, while C2-M shows better fit 
for V c > 140km/s. This difference between C2-H and C2- 
M is purely due to the fluctuation in the initial model. 
The high-wave number initial fluctuations we added for 
re-simulation are not exactly the same for these runs. 



As we have seen in figure 2, the number count for this 
run drops off from the power-law line also around V c = 
140km/s. Thus, the actual number of galaxies observed 
in the Virgo cluster goes even lower than the simulation 
which gives the number lower than the power-law. For 
reference, we drew the power-law curve in figure 5. The 
curve implies that, ideally, for V c ~ lOOkm/s, the difference 
is already significant, and for 50km/s, there is more than 
a factor of 10 difference. 

3.3. Distribution of subhalos in M — V c plane 

In figure 6 the mass of each subhalo is plotted against its 
circular velocity for runs Gl-H and Cl-H. Figure 6 shows 
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Fig. 3. Top: Cumulative number of subhalos as a function of M. Bottom: Cumulative number of M normalized by M . Left and 
right of both are for Gl halos and CI halos. Meanings of the curves are the same as in figure 2. 



that the correlation between M and V c is fairly tight, and 
well fitted by M oc V c 3 , for both runs Gl and CI. This 
tight relation, however, becomes somewhat loose for the 
velocity smaller than 20 km/s in the case of the galaxy 
scale Gl and than 100 km/s in the case of the cluster 
scale CI. We show the la unbiased variances in figure 
6, from which we can sort of sec that the distribution of 
circular velocity for a given mass range is wider for lower 
mass. 

To see this tendency more clearly, in figure 7 we plot 
the unbiased variance of the circular velocity as the func- 
tion of halo size. Here, we used the number of particles in 
subhalos as the indicator of the size of the halo. We can 
see that the variance shows essentially the same behav- 



ior for three runs with different mass resolutions, at least 
for cluster-scale runs Cl-X. This behavior indicates that 
the increase of the variance in low-mass halos is not the 
physical reality but the numerical artifact. The real vari- 
ance is probably around 15%, independent of the mass of 
subhalos. The galaxy-scale runs show a similar tendency. 

The number of particles below which the variance be- 
comes large is around 200, the same as the number of 
particles below which the number count becomes unreli- 
able. Thus, these deviations are probably driven by the 
same mechanism. 
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Fig. 4. Cumulative number of subhalos as a function of n p , 
with the fitting function of eq. (5). Bold curves are in Gl 
halos and thin curves are in CI halos. Line types indicate the 
resolution of runs, the same as figure 2 

4. Summary and discussions 

4-1. Distribution of subhalos 

We carried out large scale iV-body simulations of for- 
mation of CDM halos. We simulated the formation of 
both a galactic size halo and a cluster size halo. We de- 
termined the masses and circular velocities of subhalos for 
both cases. We found that the cumulative distribution of 
subhalos, both as a function of the circular velocity and as 
that of the mass, is affected by mass resolution of the sim- 



Fig. 5. Cumulative number of subhalos as a function of V c 
in runs C2-X. Solid, dashed and dotted curves indicate runs 
H, M, and L, as in figure 2. Smooth solid curve shows V c -3 . 
Crosses are the distribution of galaxies in Virgo Cluster, used 
in figure 2 of M99. 

ulation. In all runs, the number of subhalos n with more 
than 200 particles is well expressed by a single power-law 
of n oc V~ 3 . For the subhalos with less than 200 parti- 
cles, the number of subhalos becomes smaller than this 
power-law line. 

We investigated the distribution of subhalos in the plane 
of their masses and circular velocities, and found that the 
there is a tight correlation with M oc V 3 . The tightness of 
the correlation again depends on the number of particles 
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Fig. 7. Top: The ratio between the dispersion(cr) and the average (V c ) of V c as a function of n p . Vertical error bars are derived 
from dispersion of cr and Vc- Bottom: Average of V c as a function of M. Upper horizontal axis shows the corresponding n p for high 
resolution run. Vertical error bars indicate the dispersion of V Ct the same data as plotted in figure 6. Left and right are for Gl and 
CI runs, respectively. The meanings of curves are the the same as that in figure 2. 



in subhalos, and for subhalos with less than 200 particles 
the correlation becomes weaker. 

4--2. "Dwarf galaxy problem" in clusters of galaxies 

From runs Gl-X, the number of subhalos with n p > 200 
are about 30 for run Gl-L, 200 for run Gl-H. So roughly- 
speaking, our result implies that the cumulative distribu- 
tion of subhalos obtained by iV-body simulation is reliable 
only for N/{2 x 10 4 ) most massive halos, where N is the 
number of particles in the virial radius of the parent halo. 



Thus, results obtained with, for example, 1M particles is 
okay for only the first 50 subhalos. This result means 
that the good agreement of the observation and simula- 
tion result obtained with around 1M particles is actually 
a numerical artifact, and there is a serious discrepancy 
between the observed number of galaxies in clusters like 
the Virgo cluster and the number of subhalos obtained by 
iV-body simulation. Ideally the number of subhalos with, 
say, circular velocity of lOOkm/s is factor two or three 
times bigger than the observed count, and for 50km/s this 
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factor is more than 10. 

Most "solutions" so far proposed for the dwarf galaxy 
problem were aimed at reducing the number of dwarfs 
with V c ~ 20km/s or less. (Colin ct al. 2000, Dave ct al. 
2001, Stoehr et al. 2002, K04) However, most of them do 
not work for halos with mass more than 10 11 M Q . 

One serious theoretical problem with the subhalos with 
mass lO n M0 is that when they were first formed they 
were generally much more massive. (De Lucia et al. 2004, 
K04) Most of their mass have been stripped out through 
the tidal interaction with the parent halo. In other words, 
they were initially almost as massive as the halo of our 
galaxy. Thus, a "solution" which can account for the dis- 
crepancy between the observed number of galaxies in clus- 
ters and the subhalos in a CDM halo need to be able to 
reduce the number of halos with initial virial mass as large 
as 10 12 M Q by a factor of two or so. 

Acknowledgment We thank Toshiyuki Fukushigc 
and Hideki Yahagi for helpful discussions. The simula- 
tion was partly done using IBM pSeries 690 of the Data 
Reservoir Project of the University of Tokyo. This re- 
search is partially supported by the Special Coordination 
Fund for Promoting Science and Technology (GRAPE-DR 
project), and a grant for the 21st Century COE program 
of the University of Tokyo, "Quantum Extreme Systems 
and Their Symmetries" Ministry of Education, Culture, 
Sports, Science and Technology, Japan. 

References 

Bardeen, J. M., Bond, J. R., Kaiser, N., & Szalay, A. S. 1986, 

ApJ, 304, 15 
Barnes, J. & Hut, P. 1986, Nature, 324, 446 
Bertschinger, E. 2001, ApJS, 137, 1 

Colin, P., Avila-Reese, V., & Valenzuela, O. 2000, ApJ, 542, 
622 

Dave, R., Spergel, D. N., Steinhardt, P. J., & Wandelt, B. D. 

2001, ApJ, 547, 574 
Davis, M., Efstathiou, G., Frenk, C. S., & White, S. D. M. 

1985, ApJ, 292, 371 
De Lucia, G, Kauffmann, G, Springel, V., White, S. D. M., 

Lanzoni, B., Stoehr, F., Tormen, G, & Yoshida, N. 2004, 

MNRAS, 348, 333 
Diemand, J., Moore, B., & Stadel, J. 2004a, MNRAS, 352, 535 
Diemand, J., Moore, B., Stadel, J., & Kazantzidis, S. 2004b, 

MNRAS, 348, 977 
Fukushige, T., Makino, J., & Kawai, A. 2005, PASJ, 57, 1009 
Gao, L., White, S. D. M., Jenkins, A., Stoehr, F., & Springel, 

V. 2004, MNRAS, 355, 819 
Ghigna, S., Moore, B., Governato, F., Lake, G, Quinn, T., & 

Stadel, J. 2000, ApJ, 544, 616 
Kawai, A., Makino, J., & Ebisuzaki, T. 2004, ApJS, 151, 13 
Klypin, A., Gottlober, S., Kravtsov, A. V., & Khokhlov, A. M. 

1999a, ApJ, 516, 530 
Klypin, A., Kravtsov, A. V., Valenzuela, O., & Prada, F. 

' 1999b, ApJ, 522, 82 
Kravtsov, A. V., Gnedin, O. Y., & Klypin, A. A. 2004, ApJ, 

609, 482 
Makino, J. 2004, PASJ, 56, 521 

Makino, J., Fukushige, T., Koga, M., & Namura, K. 2003, 
PASJ, 55, 1163 



Moore, B., Ghigna, S., Governato, F., Lake, G, Quinn, T., 

Stadel, J., & Tozzi, P. 1999, ApJL, 524, L19 
Moore, B., Katz, N., & Lake, G. 1996, ApJ, 457, 455 
Navarro, J. F., Frenk, C. S., & White, S. D. M. 1996, ApJ, 
462, 563 

Pfitzer, D., Salmon, J., & Sterling, T. 1997, CACR Technical 
Reports, 150 

Reed, D., Governato, F., Quinn, T., Gardner, J., Stadel, J., & 
Lake, G. 2005, MNRAS, 359, 1537 

Spergel, D. N. & Steinhardt, P. J. 2000, Physical Review 
Letters, 84, 3760 

Springel, V., White, S. D. M., Jenkins, A., Frenk, C. S., 
Yoshida, N., Gao, L., Navarro, J., Thacker, R., Croton, D., 
Helly, J., Peacock, J. A., Cole, S., Thomas, P., Couchman, 
R, Evrard, A., Colberg, J., & Pearce, F. 2005, Nature, 
435, 629 

Stoehr, F., White, S. D. M., Tormen, G, & Springel, V. 2002, 

MNRAS, 335, L84 
White, S. D. M. & Rees, M. J. 1978, MNRAS, 183, 341 



